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INTRODUCTION 



Binary alloys involving equal proportions of a noble metal 
Cu and a divalent metal Zn have a stable low temperature 
/3-phase which sits immediately to the right of the pure 
face-centered cubic Cu phase in the alloy phase diagram 
P, 0] ■ This phase called f3 brass has a B2 cubic structure 
with two atoms per unit cell. At high temperatures the 
alloy forms a disordered body-centered cubic structure. 
The disorder-order transition takes place around 730K. 
The alloy satisfies the Hume-Rothery rules and has 
the same ratio of valence electrons to atoms. Jona and 
Marcus Q have shown from a Density Functional The- 
ory (DFT) based approach that within the Local Den- 
sity Approximation (LDA), it is the body-centered based 
B2 which is the stable ground state. They also showed 
that if we include the Gradient Corrections (GGA) then 
we get a tetragonal ground state lower in energy by 0.1 
mRy/atom. This is in contradiction with the latest ex- 
perimental data. The alloying of face-centered cubic Cu 
with an equal amount of Zn leads to a body-centered 
stable phase. Zn has only one more electron than Cu. 
This is an interesting phenomenon. CuZn alloys also 
have anomalously high elastic anisotropy. This makes 
the theoretical study of CuZn an interesting exercise for 
a proposed theoretical technique. The phase diagram of 
CuZn is shown in Fig. ^ 

One of the earliest first-principles Density functional 
based study of the electronic properties of CuZn as by 
Bansil 0- The authors had studied the complex bands 
of cv-phase of CuZn using the Korringa-Kohn-Rostocker 
(KKR) method coupled with the coherent potential ap- 
proximation (CPA) to take care of disorder. They com- 
mented on the effects of charge transfer and lattice con- 
stants on the electronic structure. They found the elec- 
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FIG. 1: Experimental phase diagram of CuZn alloy 



tronic distribution of this alloy to be of a split band kind 
with the centers of the Cu and Zn d-bands well sepa- 
rated from each other. Their Zn d-bands showed hardly 
any dispersion and were shown only schematically in their 
figures. In a later work Rowlands |5j generalized the CPA 
to a non-local version (NL-CPA) and studied the effects 
of short-ranged ordering in CuZn. Their technique was 
based on an idea of renormalization in reciprocal space 
suggested by Jarrel and Krishnamurthy 
The order-disorder transition in CuZn is a classical ex- 
ample of a true second order transition. Several very 
early investigations on this alloy has been reviewed by 
Nix and Shockley Q and Guttman |j. These investi- 
gations, of course, were rather crude, as sophisticated 
approaches to deal with disordered alloys had not really 
been developed at that time. However, it was recognized 
that a knowledge of the short-ranged order correlations 
above the critical temperature should be of considerable 
interest. Early neutron scattering experiments were car- 
ried out on /3-brass by Walker and Keating @. The 
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Warren-Cowley short-ranged order parameter, defined by 
a(R) — 1 — Pab (R) I x where x is the concentration of A, 
Pab {P) is the probability of finding an A atom at a dis- 
tance of R from a B atom, was directly obtained from 
the diffuse scattering cross-section : 

% = x(l-x)(b A -b B ) 2 Ya(R)f(K)exp(iK-R) 

where bA, bs are the scattering lengths of A and B atoms, 
and f(K) = exp(—C\K\ 2 ) is the attenuation factor aris- 
ing from thermal vibrations and static strains. The ex- 
perimental data for the short-ranged order parameter as 
a function of temperature are thus available to us. A 
Ising-like model using pair interactions was studied by 
Walker and Chipman [lfj and the short-ranged order 
was theoretically obtained. However, the pair interac- 
tions were simply fitted to the experimental values of the 
transition temperature Tq and in that sense it was an 
empirical theory. The experimental estimate of the near- 
est neighbour Warren-Cowley parameter was found to be 
varying between -0.171 to -0.182 at around 750K. 

In a later work using the much more sophisticated locally 
self-consistent Green function (LSGF) approach based on 
the tight- binding linear muffin-tin orbital (TB-LMTO) 
technique Abrikosov et.al. 0] studied CuZn alloys. The 
authors argued that earlier studies of the mixing en- 
thalpies of CuZn using the standard coherent poten- 
tial approximation approaches [l2T|-|l6l| showed signifi- 
cant discrepancy with experiment. Part of the discrep- 
ancy was assumed to be from neglect of charge transfer 
effects and part from the effects of short ranged order- 
ing (SRO). The main thrust of this technique, which was 
based on an earlier idea of a locally self-consistent mul- 
tiple scattering (LSMS) by Wang et.al. was to go 
beyond the CPA and include the effects of the immedi- 
ate environment of an atom in the solid. The LSMS gave 
an excellent theoretical estimate of the ordering energy 
in CuZn : 3.37 MRy/atom as compared to the exper- 
imental value of 3.5 MRy/atom. The LSGF approach 
correctly predicted ordering tendency in CuZn on lower- 
ing temperature and combing with a Cluster Variation- 
Connolly Williams (CVM-CW) obtained a value of the 
nearest neighbour Warren-Cowley SRO parameter a = 
-0.15. Subsequently Bruno et.al. 01 proposed a mod- 
ification of the CPA taking into account the local field 
and showed that charge transfer effects can be taken into 
account accurately as compared to the O(N) methods 
just described. They applied their approach to the CuZn 
alloys. 

One of the earlier works on the optical property of CuZn 
alloy was the determination of the temperature variation 
of spectral reflectivity by Muldawer [ijj . The author at- 
tempted to explain the color of the disordered /3-brass 
CuZn alloy via the internal photoelectric effect |20j . Al- 



though, the experimental data also contained the contri- 
bution from plasma oscillations, the author claimed that 
the spectral reflectivity helps to explain the band picture 
of the alloys as a function of the inter-atomic spacin g. I n 
order to explain the optical properties, Amar et.al. |21| - 
[23| studied the band structure of CuZn using the KKR 
method. However, they had used the virtual crystal ap- 
proximation, replacing the random potential seen by the 
electrons by an averaged one. This is now known to be 
particularly inaccurate for split band alloys. 

The above discussion was necessary to bring into focus 
the following points : in the study of alloys like CuZn it 
would be interesting to address the effects of charge trans- 
fer and short-ranged ordering. In this communication we 
shall address exactly these two points. We shall propose 
the use of the augmented space recursion (ASR) cou- 
pled with the tight-binding linear muffin-tin orbitals ba- 
sis (TB-LMTO) H3] to study the effects of short-ranged 
ordering on both the electronic structure and the opti- 
cal properties of /3-CuZn alloy at 50-50 composition. We 
should like to stress here that the TB-LMTO- ASR ad- 
dresses precisely these effects with accuracy : the Den- 
sity functional self-consistent TB-LMTO takes care of 
the charge transfer, while the local environmental effects 
which are essential for the description of SRO are dealt 
with by the ASR. The TB-LMTO- ASR and its advan- 
tages has been extensively discussed earlier by Mooker- 
jee H3 and in a series of articles [2I El EtI El El l30| . 
We would like to refer the interested readers to these for 
details. 



SPECTRAL FUNCTIONS, COMPLEX BANDS 
AND DENSITY OF STATES FOR 50-50 CUZN 

In this section we shall introduce the salient features of 
the ASR which will be required by us in our subsequent 
discussions. 

We shall start from a first principle TB-LMTO set of 
orbitals j^, H3| in the most-localized a representation. 
This is necessary, because the subsequent recursion re- 
quires a sparse representation of the Hamiltonian. In 
this representation, the second order alloy Hamiltonian 
is given by, 

H (2) = E„ + h - hoh 

where, 

h =]T (C R - E vR ) V R + Y,Ys A R 2 S ««' A R' 2 T ««' 

R R R' 

°=J2 orVr w 

R 
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C_r, E„jj , Afl and or are diagonal matrices in angu- 
lar momentum space, and is the structure matrix. 
These matrices are of rank L max . Vr = \R)(R\ and 
Trr> = \R)(R'\ are projection and transfer operators 
in the Hilbert space H spanned by the tight-binding ba- 
sis {|-R)}- Here, R refers to the position of atoms in the 
solid and L is a composite label {£, m, m s } for the angu- 
lar momentum quantum numbers. C, A and o are the 
potential parameters of the atoms which sit on the lat- 
tice sites ; o _1 has dimension of energy and E„'s are the 
reference energies about which the muffin-tin orbitals are 
linearized. 

For a disordered binary alloy we may write : 

Crl = Ct n R + Cf (1 - n R ) 
A^ 2 = (Af) 1/2 n fl+ (Af) 1/2 (l-^) 
orl = of n R + of (l-n R ) (2) 

Here {ur} are the random site-occupation variables 
which take values 1 and depending upon whether the 
muffin-tin labelled by R is occupied by an A or a B- 
type of atom. The atom sitting at {R} can either be of 
type A {riR = 1) with probability x or B (jir = 0) with 
probability y, where x and y are the concentrations of 
the components A and B in the binary alloy. The aug- 
mented space formalism (ASF) now introduces the space 
of configurations of the set of binary random variables 
{tir} : <&. In the absence of short-ranged order, each 
random variable ur has associated with it an operator 
Mr whose spectral density is its probability density : 

p(n R ) = x8(n R - 1) + yS(n R ) 

= -- km Im(T* | ((ur + iS)I - Mr)- 1 \ U)- 

7T S— >0 

(3) 

Mfl is an operator whose eigenvalues 1,0 correspond to 
the observed values of ur and whose corresponding eigen- 
vectors {|1_r), \0r)} span a configuration space 4>r of rank 
2. We may change the basis to {| | Ir)} 

\U) = {V^\o R ) + Vy\iR)} 
llfl) = {Vy|Oft)-Vi|ifi>} 

and in the new basis the operator Mr corresponding to 
n R is : 

Mr = xV R + yV R + ^yTl [ 
= xl+b(x)V l R + f(x)T^ 

where v\ = \U)(Ul K = \Ir)(Ir\ and = 
Wr)(Ir\ + Ufl)(Tfll are the projection and transfer oper- 
ators in the configuration space 4>R spanned by the two 
basis vectors. b(x) = (y — x) and f(x) = y/xy. 



The full configuration space $ = Y[ R 4>R 1S then spanned 
by vectors of the form | Tt-Lt-L ■ ■ ■)■ These configurations 
may be labelled by the sequence of sites {C} at which 
we have a |. For example, for the state just quoted {€} 
= | {3, 5,...}). This sequence is called the cardinality 
sequence. If we define the configuration | ff . . . f . . .) as 
the reference configuration, then the cardinality sequence 
of the reference configuration is the null sequence {0}. 

In the full augmented space the operator corresponding 
to Ur is : 

M fl = iI®I... + b(x) V R ®1... + f(x) Tl [ <g> J . . . 

(4) 

The augmented space theorem [2(| states that 

« A({n R }) » = < m\A({M R })\{(t)} >, (5) 

where, 

A({M R }) = J...jA({X R })l[dP(X R ). 

P(Aij) is the spectral density of the self-adjoint operator 
Mr. Applying JSJ) to the Green function we get : 

«G(k, 2 )»= (k® {0}|(zI-H (2) ) _1 |k® W)- (6) 

where G and H^ 2 ^ are operators which are matrices in 
angular momentum space, and the augmented reciprocal 
space basis |k, L ® {0}) has the form 

(1/VN) exp(-ik ■ R)\R 7 L® {0}). 

R 

The augmented space Hamiltonian H' 2 ) is constructed 
from the TB-LMTO Hamiltonian H' 2 ' by replacing each 
random variable ur by the operators Mr. It is an oper- 
ator in the augmented space \& = TL ® <&. The ASF maps 
a disordered Hamiltonian described in a Hilbert space 
7i onto an ordered Hamiltonian in an enlarged space \f, 
where the space ^ is constructed as the outer product 
of the space Ti and configuration space $ of the random 
variables of the disordered Hamiltonian. The configu- 
ration space 4> is of rank 2 N if there are N muffin-tin 
spheres in the system. Another way of looking at H^ 2 ^ is 
to note that it is the collection of all possible Hamiltoni- 
ans for all possible configurations of the system. 

Combining (2) and (4) we get for any of the operators V 
diagonal in real and angular momentum spaces : 
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A + B 



(7) 



r 



where for any of the diagonal (in real and angular mo- 
mentum spaces) operator V : 



A(V)= (xV£ + yVf)6 
B(V)= b(x)(Vf-V L B )S 
F(V)= f(x){V£-V*)5 



LL' 



LL' 



LL' 



(8) 



In case there is no off-diagonal disorder due to local lat- 
tice distortion because of size mismatch : 



R R' 

This equation is now exactly in the form in which the 
recursion method |33l | may be applied. At this point we 
note that the above expression for the averaged GllO*., z) 
is exact. 

The recursion method addresses inversions of infinite ma- 
trices. Once a sparse representation of an operator in 



Hilbcrt space, H' 2 '. is known in a countable basis, the 
recursion method obtains an alternative basis in which 
the operator becomes tridiagonal. This basis and the 
representations of the operator in it are found recursively 
through a three-term recurrence relation : 

\u n+1 } = H^K} ~ a n (k)\u n } - 0l(k)\u n ^}. (9) 

with the initial choice \ui} = \kL) g) |{0}) and 0\ = 1. 
The recursion coefficients a n and (3 n are obtained by im- 
posing the ortho-normalizability condition of the new ba- 
sis set as : 



a„(k) = 



{n|H( 2 )|n} 



{n-l|H( 2 )|n} 



{n\n} ' ^"" iv ' {n\n} 
and also {mlH^ 2 -* |n} = for m ^ n, n ± 1 

To obtain the spectral function we first write the configu- 
ration averaged L-projected Green functions as continued 
fractions : 



J 



■C G L i(k,z) > 



z - an,(k) - 



z - a 2 L(k) 



/? 3 2 L(k) 



r 



z - ajvi(k) - r L (k, z) 



where T j, (k, z) is the asymptotic part of the continued 
fraction. The approximation involved has to do with 
the termination of this continued fraction. The coeffi- 
cients are calculated exactly up to a finite number of 
steps {a n , f3 n } for n < N and the asymptotic part of the 
continued fraction is obtained from the initial set of co- 
efficients using the idea of Beer and Pettifor terminator 
[34| . Haydock and coworkers j3f| have carried out exten- 
sive studies of the errors involved and precise estimates 
are available in the literature. Haydock [36| has shown 
that if we carry out recursion exactly up to N steps, the 
resulting continued fraction maintains the first 2N mo- 
ments of the exact result. 

It is important to note that the operators A, B,F are 



all projection operators in real space (i.e unit operators 
in k- space) and acts on an augmented space basis only 
to change the configuration part (i.e. the cardinality se- 
quence {C} ). 



A||{C}> - A 1 \\{C}), 

B\\{C}) = A 2 \\{C})S(R€{C}), 

F\\{C}) = A 3 \\{C±R}). 

The coefficients A\ — A3 can been expressed from equa- 
tion JSJ). The remaining operator S is diagonal in k- 
space and acts on an augmented space only to change 
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the configuration part : 



S||{C}}=£exp(-*k. X )||{C-x}}. 



Here x s are the near neighbour vectors. The operation 
of the effective Hamiltonian is thus entirely in the con- 
figuration space and the calculation does not involve the 
space H, at all. This is an enormous sim plif ication over 
the standard augmented space recursion 0, EE EE EE| ; 
where the entire reduced real space part as well as the 
configuration part were involved in the recursion process. 
Earlier we had to resort to symmetry reduction of this 
enormous space in order to make the recursion tractable. 
Here the rank of only the configuration space is much 
smaller and we may further reduce it by using the local 
symmetries of the configuration space, as described in 
our earlier letter E3- However, this advantage is offset 
by the fact that the effective Hamiltonian is energy de- 
pendent. This means that to obtain the Green functions 
we have to carry out the recursion for each energy point. 
This process is simplified by carrying out recursion over 
a suitably chosen set of seed energies and interpolating 
the values of the coefficients across the band. 

The self-energy which arises because of scattering by the 
random potential fluctuations is of the form : 



S L (k,z) 



/? 2 2 L(k) 



"2L(k) 



a NL (k) - T L (k, z) 



So the continued fraction can be written in the form 
l/{z - E l {\l) - V l (\<l,E)), where E L (k) = a 1L (k). 

The average spectral function -c A k {E) 3> is related to 
the averaged Green function in reciprocal space as : 

« A k (E) > = < A kL (E) », 



where 

< A kL (E) > = -- lim {Im < G LL (k,E -id) >}. 

To obtain the complex bands for the alloy we fix a value 
for k and solve for : 



z-Si(k)-E L (k,£?)=0. 



The real part of the roots will give the position of the 
bands, while the imaginary part of roots will be propor- 
tional to the lifetime. Since the alloy is random, the 
bands always have finite lifetimes and arc fuzzy. 



We have used this reciprocal space ASR to obtain the 
complex bands and spectral functions for the CuZn al- 
loy. This is shown in Figs. © - ©■ It should be noted 
that we have carried out a fully LDA self-consistent cal- 
culation using the TB-LMTO-ASR developed by us Ell 
to obtain the potential parameters, ft takes care of the 
charge transfer effects. For the Madelung energy part 
of the alloy calculation, we have chosen the approach of 
Ruban and Skriver |3sj . 

The two panels of Fig. J5J) compare the band structures 
of pure Cu and pure Zn metals in the same bec lattice 
as the 50-50 alloy. We note that the s-like bands of Cu 
and Zn stretch from -0.8 Ry, while the d-like states of 
Zn and Cu, whose degeneracies are lifted by the cubic 
symmetry of the bec lattice are more localized and reside 
in the neighbourhood of -0.6 Ry. and between -0.3 to - 
0.2 Ry. respectively. The complex bands of the solid 
clearly reflect the same band structure. However, the 
bands are slightly shifted and broadened because of the 
disorder scattering of Bloch states in the disordered alloy. 
The scattering lifetimes are maximum for the Cu (i-likc 
bands, less for the Zn <i-likc bands and minimum for the 
lower s-like bands. This is expected, since the delocalized 
s-like states straddle large volumes of the lattice space 
and are therefore less sensitive to the local configuration 
fluctuations of the substitutional alloy. 

The same is reflected in the spectral functions, shown 
here also along the T-X direction in the Brillouin zone. 
Sharp peaks stretching from -0.8 Ry, groups of wider 
peaks around -0.6 Ry. with less dispersion, characteris- 
tic of the more localized e?-like states and groups of much 
wider peaks straddling -0.3-0.2 Ry. also with less disper- 
sion. The spectral functions play an important role in 
response functions related to photoemmission and opti- 
cal conductivity |58j . Our complex bands agree remark- 
ably well with the Fig. 3 of Bansil and Ehrenreich p|. 
These authors of course did not show the dispersion of 
the Zn d-bands, but as in their work, the Cu bands show 
greater disorder induced lifetimes than the lower energy 
Zn d-bands. 

We may use the generalized tetrahedron method to pass 
from the reciprocal space spectral functions to the real 
space density of states . Alternatively, we may also 
carry out real-space ASR to obtain the density of states 
directly. 

Fig. |J3J shows the densities of states for the pure Zn 
(solid lines) and Cu (dashed lines) in the same bec lat- 
tice as the alloy and compares this with the ordered B2 
and disordered bec 50-50 CuZn alloy. We first note that 
in the ordered B2 alloy there is a considerable narrow- 
ing of the Zn well as the Cu d-like bands. The feature 
around -0.3 Ry. is suppressed in the ordered alloy. In 
the disordered alloy on the other hand, although disor- 
der scattering introduces life-time effects which washes 
out the sharp structures in the ordered systems, the re- 
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FIG. 2: (left) Bands for pure Cu and Zn in bcc lattices with the same lattice parameter as the 50-50 CuZn alloy. The dashed 
line are for Cu and the full lines for Zn. (right) Complex bands for the 50-50 CuZn alloy. 



semblance to the pure metals is evident. As seen in the 
complex bands, the life-time effects in the Cu d-like part 
is prominent. If we interpret the left most figure as that 
due to completely segregated Cu-Zn and the middle one 
as the completely ordered one, then the disordered alloy 
lies between the two. In the next section, introducing 
short-ranged ordering effects on top of the fully disor- 
dered alloy, we shall study how to bridge between the 
two states. 




-0.6 -0.4 
Energy (Ryd.) 



FIG. 3: Spectral functions for the CuZn alloy for k-vectors 
along the r to X direction in the Brillouin zone. 



SHORT-RANGED ORDERING IN THE ALLOYS 

Attempts at developing generalizations of the coher- 
ent potential approximation (CPA) to include effects of 
short-ranged order (SRO) have been many, spread over 
the last several decades. The CPA being a single site 
mean-field approximation cannot take into account SRO 
as it spans at least a nearest neighbour cluster on the 
lattice. The early attempts at cluster-generalizations of 
the CPA were beset with difficulties of violation of the 
analytic properties of the approximated configuration av- 
eraged Green function. Tsukada's idea of introducing 
a super-cell the size of the cluster in an effective medium 
suffered from the problem of broken translational symme- 
try within the cluster even when the disorder was homo- 
geneous. This same problem beset the CCPA's proposed 
by Kumar et.al. [^3 based on the augmented space the- 
orem introduced by Mooke rjee |2r3 | . The embedded clus- 
ter method of Gonis et.al. |42( immersed a cluster in a 
CPA medium which lacked the full self-consistency with 
it. The first, translationally symmetric cluster approx- 
imations which preserved the analytic properties of the 
approximate Green functions were all based on the aug- 
mented space theorem of Mookerjee [2(|. They included 
the travell ing cluster approximation (TCA ) o f Kaplan 
and Gray |2jj and Mills and Ratnavaraksa |44j and the 
CCPA proposed by Razee et.al. j^. The problem with 
these approaches was that they became intractable as 
the size of the cluster considered was increased much be- 
yond two sites. Mookerjee and Prasad ^(| generalized 
the augmented space theorem to include correlated dis- 
order. However, since they then went on to apply it in 
the CCPA approximation, they could not go beyond the 
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FIG. 4: (left) Density of states of pure Zn (solid lines) and Cu (dashed lines)in the same bcc lattice as the 50-50 CuZn alloy, 
(centre) Density of states for ordered B2 50-50 CuZn alloy, (right) Density of states for the disordered bcc 50-50 CuZn alloy. 




FIG. 5: (Color Online) (left) The blue atom represents the site labelled 0, where the local density of states will be calculated. 
The red atoms are its eight nearest neighbours on the bcc lattice. The white atoms are more distant neighbours (right) 
Schematic diagram showing lattice translational symmetry of the nearest neighbour cluster SRO approach. 



two-site cluster and that too to model systems alone. The 
breakthrough came with the augmented s pac e recursion 
(ASR) approach proposed by Saha et.al. [i^-jisl]. The 
method was a departure from the mean-field approaches 
which always began by embedding a cluster in an effective 
medium which was then obtained self-consistcntly. Here 
the Green function was expanded in a continued fraction 
whose asymptotic part was obtained from its initial steps 
through an ingenious termination procedure ■ In this 
method the effect at a site of quite a large environment 
around it could be taken into account depending how far 
one went down the continued fraction before termination. 
The technique was made fully LDA-self-consistent within 
the tight-binding linear muffin-tin orbitals (TB-LMTO) 
approach |4flj and several applications have been carried 
out to include short-ranged order in different alloy sys- 
tems |5{J . Recently Leath and co-workers have developed 
an itinerant CPA (ICPA) based on the augmented space 
theorem j^jj, which also maintains both analyticity and 
translational symmetry and takes into account effect of 
the nearest neighbour environment of a site in an alloy. 
The technique has been very successfully applied to the 
phonon problem in alloys where there were large force 
constant disorders. This method matches well with the 
ASR approach to the same alloys and there is now a 
concerted effort to apply it to electronic problems based 



on both the TB-KKR and the TB-LMTO methods. A 
very different and rather striking approach has been de- 
veloped by Rowlands et.al. (the non-local CPA or 
NL-CPA) using the idea of coarse graining in reciprocal 
space originally proposed by Jarell and Krishnamurthy 
The NL-CPA with SRO has been applied earlier by 
Rowlands et.al. [j| and is on the verge of being made fully 
DFT self-consistent within the KKR. The authors report 
an unpublished report on it |54j . fn this communication 
we report a fully DFT self-consistent ASR based on the 
TB-LMTO with SRO incorporated. We have applied it 
to the case of 50-50 CuZn alloys, so as to have a compar- 
ison with earlier attempts using different techniques. 



THE GENERALIZED AUGMENTED SPACE 
THEOREM 

The generalized augmented space theorem has been de- 
scribed in detail by Mookerjee and Prasad Let us 
briefly introduce those essential ideas which are necessary 
to make this communicated reasonably self-contained. 

For a substitutionally binary disordered alloy A x B y on a 
lattice we can introduce a set of random occupation vari- 
ables {rift} associated with the lattice sites labelled by 
i?, which take the values or 1 depending upon whether 
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the site R is occupied by a A or a B type of atom. The 
Hamiltonian and hence the Green function arc both func- 
tions of this set of random variables. The configuration 
averaged Green operator is then : 



«GO,{n R })>= /••■/ Y[dm 



■G(z,{n R })P({n R }) 



To start with, let us assume that short-ranged order ex- 
tends up to nearest neighbours only. Let us take for an 
example the nearest neighbour cluster of nine atoms on 
a body centered cubic lattice (see Fig. [SJ) centered on the 
site labelled by i = 0. Let us concentrate on the central 
atom (labelled 0). The occupation variables associated 
with its eight neighbours arc correlated with n\ , but not 
with one another. Further none of the other occupation 
variables associated with more distant sites are correlated 
with n\. If we label any other site k as 0, that is carry 
out a lattice translation from the site i to k, its environ- 
ment is identical to the earlier situation. Therefore the 
subsequent procedure is translationally symmetric, as it 
should be, if the SRO is itself homogeneous. We may 
then write : 



While for j=l,2,. . . 8 : 

P(n,j\no = 1) = (x + ay)5(rij — 1) + (1 — a)yS(rij) 
P(n,j\no = 0) = (1 — a)xS(rij — 1) + (y + ax)6(rij) 

(10) 

In terms of the concentrations x and y and SRO param- 
eter a these operators have representations : 



M k = 
M) = 
M° = 



xy y 
x + ay 



0, or k > 8 



V(l -a)y(x + ay) 
y/ '(1 - a)y{x + ay) (1 - a)y 



(1 — a)x — a)x(y + ax) 

3 \ — a)x{y + ax) y + ax 

3 = 1,2. ..8 (11) 



P 1 



The projection operators at the site labelled are 

x y/xy 
/xy y 

V -\jxy 

-Jxy x 



(12) 



P(n ,n 2 ,...n k ...)=P(n ) ] [ P(n>i) ] [ P(n k ) 



k>8 



The generalized Augmented Space Theorem then asso- 
ciates with the random variables {n k } corresponding op- 
erators {Mfc} in their configuration space. The construc- 
tion of the representations of these operators has been 
discussed in detail in the paper by Mookerjee and Prasad 
j4t| . Here we shall quote only the relevant results neces- 
sary to proceed further. 

We shall characterize the SRO by a Warren-Cowley pa- 
rameter a. In terms of this the probability densities are 
given by : for k=0 and k>8. 

P(n k ) = x5(n k - I) + y5(n k ), x + y = l 



In the full augmented space, which is the product space 
of all the individual configuration space of the different 
occupation variables, the operators which replace the oc- 
cupation variables are : 



Mo 

M, 

M fc 



Mi g) I <g) I 1 



I 



fc=0 

I ® i" ® . . . Mfc g> I ( 



j=l,2,... 8 

k > 8 



(13) 



In operator form the operators labelled by and k > 8 
have the same form as in Equation The operators 
labelled by j have the form : 



Mj = xl®l... + b{x)V}®!...+ b'{x)V^®T... + f"(x) T U ® 1. . . + b"{x)V^®V\®... 
+ f'(x) T^®1... + d"(x) Vl ®Tj l ... + d{x) T n <g> V\ . . . + f"'(x) T U g> T? 1 . . . 

(14) 

where, 

b'(x)=a(y-x) b"(x) = -2a(y - x) F'(x) = (xfa + yfi 2 ) F"(x)=ay/xy 
F"'{x) = 09i - /3 2 )V^y; D(x) = -2a^xy; D"(x) = (y - x){pi - ft) 
Pi = a)y(x + ay) ; ft = (1 - a)x(y + ax) 



9 



We now follow the augmented space theorem and replace 
all the occupation variables {ur} by their corresponding 
operators. 

It is easy to check that in the absence of short-ranged or- 
der the functions b'{x), b"{x), f"(x), f"'(x), d(x) and 
d'(x) all vanish, (3\ = (3 2 = y/xy and f'(x) = f{x). We 
also note that the choice of the central site labelled '0' 
is immaterial. If we translate this site to any other and 
apply the lattice translation to all the sites, hence the 
Hamiltonian in the full augmented space, remains un- 
changed. This formulation of short ranged order also pos- 
sesses lattice translational symmetry, provided the short- 
ranged order is homogeneous in space. 



tendency, both the split bands narrow and lose structure. 
The feature around -0.3 Ry. disappears. This band nar- 
rowing and suppression of the feature around -0.3 Ry. is 
clearly seen in the ordered B2 alloy shown in Fig. Q). 

Our analysis is closely similar to that of Rowlands et.al. 
0. However, there are quantitative differences as both 
the basic technique (KKR and TB-LMTO) and the ap- 
proximations are quantitatively different. Finally, in Fig. 
IHlwe show the band energy as a function of the nearest 
neighbour Warren-Cowley parameter. The minimum oc- 
curs at the ordering end, as expected. Experimentally 
the alloy does show a tendency to order at lower temper- 
atures. 



EFFECT OF SRO ON THE DENSITY OF STATES 

We have carried out the TB-LMTO-ASR calculations on 
CuZn with a lattice constant of 2.85 A. The Cu and Zn 
potentials are self-consistently obtained via the LDA self- 
consistency loop. All reciprocal space integrals are done 
by using the generalized tetrahedron inte gra tion for dis- 
ordered systems introduced by us earlier [5JJ. 

To discuss the effect of SRO, leading, on one hand, to 
ordering (a < 0) and segregation on the other (a > 0), 
let us first look at Figs. |21 and 0] The complex band 
structure shown in Fig. [21 shows that the system is a 
split band alloy. The positions of the d-bands of Cu and 
Zn are well separated in energy. This implies that the 
"electrons travel more easily between Cu or between Zn 
sites than between unlike ones" So when the alloy 
orders and unlike sites sit next to each other, the overlap 
integral between the like sites decrease. This leads to a 
narrowing of the bands associated with Cu and Zn. A 
comparison between the leftmost and central panels of 
Fig. 0| shows that the bandwidths in the latter are al- 
most half of the former. This should be the main effect 
of ordering setting in. On the other hand, when the alloy 
is completely disordered, the bands gets widened by dis- 
order scattering and the sharp structures in the density 
of states are smoothened. 

Fig © shows the density of states with increasing posi- 
tive a. Positive a indicates a clustering tendency. Com- 
paring with Fig. we note that as clustering tendency 
increases the density of states begins to show the struc- 
tures seen in the pure metals in both the split bands. For 
large positive a there is still residual long-ranged disor- 
der. This causes smoothening of the bands with respect 
to the pure materials. For these large, positive as, we 
notice the development of the structure around -0.3Ry. 

Fig. J7J shows the density of states with increasing neg- 
ative a which indicates increasing ordering tendency. On 
the bec lattice at 50-50 composition we expect this order- 
ing to favour a B2 structure. With increasing ordering 



OPTICAL PROPERTIES OF CUZN ALLOYS 

In an earlier work |55| we had developed a methodology 
for the calculation of configuration averaged optical con- 
ductivity of a disordered alloy based on the augmented 
space formalism. Here we shall present the salient fea- 
tures required for the calculation of optical response func- 
tions in the 50-50 CuZn alloy. 

In linear response theory, at zero temperature, the gener- 
alized susceptibility of a disordered alloy is given by the 
Kubo formula : 



(j"(t)>=/ X>* v (t-t')A v (t') 
where, A u (t) is the vector potential, and 

x ^(t-t') = (i/h) e(*-f) <rer (*')]) 

is the current operator and O is the Heaviside step 
function. If the underlying lattice has cubic symmetry, 
X^ v = X The fluctuation dissipation theorem then 
relates the imaginary part of the Laplace transform of 
the generalized susceptibility to the Laplace transform of 
a correlation function : 

X "(uj) = (l/2ti) (1 - exp{-/3hu}}) S(u) 

where, 

S(u) = dt aip{i{uj + id)t} Tr (j"(t) j"(0) ) (15) 

Since the response function is independent of the direc- 
tion label /i for cubic symmetry, in the following we shall 
drop this symbol. In case of other symmetries we have 
to generalize our results for different directions. Our goal 
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Energy (Ry.) Energy (Ry.) 



FIG. 6: Density of states for 50-50 CuZn with increasing positive a which describes increasing clustering tendency 




will be, given a quantum "Hamiltonian" H to obtain the 
correlation function : Sit) = (0|j(i)j(O)|</>) 

We shall now determine the correlation directly via 
the recursion method as described by Viswanath and 
MiiUerj^. 

For a disordered binary alloy, S(t) = S[W.({nn})]. The 



augmented space theorem then states that : 



« s(t) »=< (MtWM) » 

= (0 ® {0}|j(t)j(O)|0 ® {0}) = S[U({M R })} 
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FIG. 8: The Band energy (from which the contribution of 
the core electrons have been subtracted) as a function of the 
nearest-neighbour Warren-Cowley SRO parameter 



We now expand the bra ('0(^)1 m this orthogonal basis : 

oo 

= £</*|5fe(*) 

k=0 

Continuing exactly as in the last section we get, 

D (]t) = m) ® {0}|/ o ) =< S(t) » (18) 

Taking Laplace transforms and using the three term re- 
currence, 

d (z) = (19) 

Pi 

— i>i 

z — a.\ 

z — a% — ... 

The configuration averaged structure function, which is 
the Laplace transform of the averaged correlation func- 
tion can then be obtained from the above : 



where the augmented space Hamiltonian and the cur- 
rent operators are constructed by replacing every random 
variable n R by the corresponding operator M fl . 

The recursion may now be modified step by step in the 
full augmented space : 



m)\ = (0®{0}iJw 

The time evolution of this bra is governed by the 
Schrodinger equation 



dt 



H 



(16) 



As before, we shall generate the orthogonal basis {(/fc|} 
for representation of equation (|16fl : 

(i) We begin with initial conditions : 

</_i| = ; (/o| = (0®{0}|I(O) 



< S(u>) > = lim 2 5Re d (u; + iS) (20) 

5-»0 

The imaginary part of the dielectric function is related 
to this correlation function through : 



e 2 (w) 



< S{uj) > 



The real part of the dielectric function is related to the 
imaginary part through a Kramers-Kronig relationship. 
The equations l|17|> - (|20|) will form the basis of our calcu- 
lation of the configuration averaged correlation function. 

Next we look at the expression for the current operator 
in the TB-LMTO basis : 



\xr) = \4>i 



w 



The dot refers to derivative with respect to energy. In 
this basis, the matrix elements of the current operator 
can be written as 



(ii) The new basis members are generated by a three 
term recurrence relationship : 

(f k+1 \ = (/ fe |H - (/ fe | & k - (/ fe _i| 0l k = 0,1,2... 
where, 



(MH|/ fc ) 

(fk\fk) 



(31 = 



(MA) 
(/fc-il/fc-i) 



(17) 



where, 



J RR r — e 



Vj&" 5 RR , + £ V^f h R " R > 



£ h ™" V R"£ + £ £ h ™» V *»WI h *"* 



R" 



R" R" 



(21) 
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with, 



RR' 

t/(3),M 
RR' 



b R 'W^ R ) ; ^ = <<^'K|<^) 



where x = -R — -R'- 

In a disordered alloy, the representations of the current 
operators are random : 



The technique for calculating these matrix elements has 
been described in detailed by Hobbs et.al. I57H . W e have 
also used this technique in our earlier paper [58j and we 
shall use it here as well. The readers are referred to these 
two papers for details. 

Ideally the next step would be to calculate 
Jaa^bb^ab^ba as the current terms between 
two sites when they are occupied by atom pairs AA, BB, 
AB and BA embedded in the disordered medium. A 
simpler first step would be to obtain these current terms 
from the pure A and B and from the ordered AB alloy. 
In general, the current operator can be written as : 



R R' 



(22) 



J"(0) 



^ B (0) (1 



Ik) 



and 



j aa(x) riR n R , + J RB (x) (1 - n R ) (1 - n R >) - 
j ab(x) n R (1 - nR') + J^ba(x) (1 - n R )n R > 



Where R—R' = \- Using the augmented space theorem, 
we replace the random variable n R by an operator M R 
in the expressions for J M (0) and J tt (x) we get the current 
operator in the augmented space : 



I" = Er {J'bbWVr ® I. ■ ■ + J?(0)Vr ® M R } + Zr Er< {J\ 

JZix) t rr , ®m r ®m w } 



RR' 



1 + Jf(x) t rr , <g> (m r + m r ^ 



(23) 



■#(x) 



-^b(o)]; #(o) = 
^b(x)]; = 

; j aa(x) - Jab(x) 



= [^a(0)-^ b (0)] 

^(xWiUx)] 

-^(X) + ^ B (X) 



The relevant operators are now chosen from Equations 
(0} and l|14|) and the current operator in full augmented 
space is set up from them. This augmented current op- 
erator is used to construct the starting state of the re- 
cursion as described in the last section. 

Since the real part of the dielectric function ei(to) can 
be now obtained from the imaginary part C2(w) using a 
Kramers Kronig relationship, we know the full complex 
dielectric function. All optical response functions may 
now be derived from these. If we assume the orientation 
of the crystal surface to be parallel to the optic axis, the 
reflectivity R{u>) follows directly from Fresnel's formula : 



R(w) = 



\AR-i 



Fig. [5] (top panel) shows the imaginary part of the dielec- 
tric function varying with frequency showing increasing 



segregating tendency. If we examine the density of states 
for both the ordered and disordered alloys we note that 
it is only around or just above hw ~ 0.1 Ry that the 
transitions from the d-bands of Cu begin to contribute 
to eaC^)- Below this energy, the behaviour is Drude like. 
This is clearly seen in the panels of the figure. As a in- 
creases and the alloy tends to segregate, since the weight 
of the structure in the density of states nearer to the 
Fermi-energy, increases, this contribution leads to the in- 
creasing weight of the structure near 0.1 Ry. 
Fig. [5] (bottom panel) shows the variation of with 
the nearest neighbour Warren-Cowley SRO parameter 
showing increased ordering tendency. The narrowing of 
the band with ordering and the increase of weightage to 
the structure away from the Fermi energy is clearly re- 
flected in the dielectric function. It is also obvious from 
the figures that the imaginary part of the dielectric func- 
tion is more sensitive to the variation of short-ranged 
order than the density of states. 
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Energy (Ry.) 

FIG. 9: (Top panel) Imaginary part of the dielectric function as a function of energy and its variation with the nearest neighbour 
Warren-Cowley SRO parameter, a > increasing, indicating segregating tendency. (Bottom panel) The same with a < 
increasing, indicating ordering tendency 



Color changes in CuZn brasses have been observed both 
with changing composition, hence different phase struc- 
tures, and on heating [J^. In particular as these alloys 
are coloured and the changes can be observed visually. 
A study of the changes in spectral reflectivity due to or- 
dering or segregation would be interesting. The colour 
of these alloys should be due to the same physical mech- 
anism of internal photoelectric excitations, as proposed 
for coloured metals like Cu and Au. Of course, we should 
take care of plasma effects and other phenomena that to- 
gether and in a combined way give rise to coJour. How- 
ever, it would be, as an initial study examine the effect 
of SRO on the spectral reflectivity of the alloy. This is 
sown in Fig. 

The first thing we note that for the completely disordered 



alloy the reflectivity has a maximum around A =5200 A 
which is not very different from that of Cu at low tem- 
peratures and is in the Yellow-Green region. Below this 
the reflectivity sharply drops (the reflectivity edge) and 
light of lower wavelengths is not reflected. With increas- 
ing segregation tendency this maximum drops almost lin- 
early with a to around 4950 A which is towards the green 
range. On the other hand with increasing ordering ten- 
dency the maximum moves out towards 5600A which 
moves towards the yellow from the yellow-green wave- 
lengths. We would expect then a shift from a yellow- 
green colour in the disordered to a more yellow region in 
the ordered material. We should note that these calcu- 
lations arc all done at OK. In an actual experiment the 
disorder-order transition takes place with varying tern- 
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FIG. 10: (Top panel) Reflectivity as a function of the wavelength and its variation with the nearest neighbour Warren-Cowley 
parameter a > increasing indicating segregating tendency. (Bottom panel) The same with a < indicating ordering 
tendency 



perature and the temperature effects on the electronic 
structure must be taken into account, as well as effect of 
plasma oscillations. We shall leave this for a later com- 
munication. 
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